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APRESENTAÇÃO 


Estimativa da Difusividade Térmica de Ligas Metálicas utilizando um Campo de 
Temperatura Periódico detalha o uso do método de Angstrom para estimar a difusividade 
térmica de materiais metálicos, especificamente para aços inoxidáveis comerciais AISI 
304 e AISI 316. Nesta técnica, a amostra é submetida a aquecimento periódico a uma 
determinada frequência térmica. A amplitude e a fase da onda térmica são medidas em 
determinados pontos do domínio físico, quando o campo de temperatura atinge o regime 
permanente periódico. Esses parâmetros são inseridos em um modelo matemático, para 
então determinar a difusividade térmica do material em análise. O sistema de medição 
construído para este trabalho foi projetado com base nos requisitos fundamentais do método 
de Angstrom. Este dispositivo é versátil, pois com algumas modificações também pode 
determinar a difusividade térmica de materiais não metálicos, como os polímeros. Além 
disso, o sistema de medição permite duas alternativas de configurações experimentais. 
Em uma delas, o lado da amostra é exposto ao vácuo e no outro é isolado termicamente. 
Ambas as configurações têm como objetivo minimizar as perdas transversais de calor, 
garantindo que a transferência de energia térmica seja unidimensional, no sentido axial da 
amostra. A difusividade térmica estimada em ambas as situações está de acordo com a 
literatura, com desvios absolutos inferiores a 12.0 %. 


A grande contribuição deste trabalho é o projeto e validação de um sistema de 
medição multifuncional e robusto, feito com materiais de baixo custo, para estimativa da 
difusividade térmica, baseado no método de Angstrom. Esta técnica é de grande interesse 
pela clareza de seu tratamento matemático e eficácia experimental, o que a torna uma boa 
ferramenta para identificar a difusividade térmica de novas ligas metálicas fabricadas em 
laboratórios de pesquisa. 


José Aguiar dos Santos Junior 


Professor Adjunto do Departamento de Engenharia Mecânica da Universidade 
Federal de Sergipe 
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RESUMO 


A caracterização termofísica dos materiais utilizados na Engenharia é de grande importância 
para realização de projetos nas mais diversas áreas de conhecimento onde os fenômenos 
ligados aos processos de transferência de calor exercem um papel fundamental. A 
difusividade térmica é uma propriedade termofísica importantíssima na análise de problemas 
de difusão de energia térmica. Este trabalho teve como proposta a determinação desta 
propriedade utilizando um campo de temperatura periódico. Para isto, foi construído um 
dispositivo experimental com princípio de funcionamento no método de Angstrom, o qual faz 
uso de um fluxo de calor periódico de uma fonte controlada, gerando assim, um campo de 
temperatura periódico na amostra em teste. Termopares foram instalados nas amostras para 
captar os sinais de temperatura gerados pelo fluxo de calor periódico. A amplitude e a fase 
destes sinais foram obtidas por meio de um software de análise gráfica. O termopar mais 
próximo da fonte de calor foi adotado como referência, ao passo que a razão de amplitudes 
e a defasagem, entre os sinais térmicos registrados pelos demais termopares em relação 
ao registrado por aquele termopar, foram calculadas. Estes resultados foram utilizados em 
modelos matemáticos para determinar a difusividade térmica, que pode ser identificada ou 
através da razão de amplitudes ou através da defasagem entre os perfis de temperatura. As 
amostras utilizadas neste trabalho foram de aço inox AISI 304 e aço inox AISI 316. Os valores 
de difusividade térmica identificados para estes materiais, quando foram comparados com 
valores disponíveis na literatura, obtiveram uma boa concordância, tendo em vista a faixa de 
incerteza apresentada. 

PALAVRAS-CHAVE: Difusividade Térmica. Propriedades Termofísicas. Campo de 
Temperatura Periódico. Método de Angstrom. Dispositivo Experimental. 


ABSTRACT 


Thermophysical characterization of materials used in engineering is very important for 
realization of projects in the most diverse areas of knowledge where the phenomena related 
to the process of heat transfer play an important role. Thermal diffusivity is a very important 
thermal property on the analysis of problems of diffusion of thermal energy. This work 
proposes the determination of this property using a periodic temperature field. For this, an 
experimental device was built with principle of operation in Angstrom's method, which makes 
use of a periodic heat flow from a controlled source, thereby generating a periodic temperature 
field in the test sample. Thermocouples were installed on the samples for capture of signals 
generated by the periodic heat flow. Amplitude and phase of these signals were obtained 
by means of graphic analysis software. The thermocouple closest to the heat source was 
adopted as reference, and the ratio and phase lag, between the thermal signals registered 
by other thermocouples in relation to that registered by that thermocouple, were calculated. 
These results were utilized in mathematical models to determine the thermal diffusivity, whose 
identification can be performed either through the amplitude ratio or through the phase lag 
between the temperature profiles. Samples utilized in this work were stainless steel AlSI 
304 and stainless steel AISI 316. Thermal diffusivity estimative, when compared whit values 
available in literature obtained a good agreement, considering the range of uncertainty 
presented. 

KEYWORDS: Thermal Diffusivity. Thermophysical Properties. Periodic Temperature Field. 
Angstrom's Method. Experimental Device. 
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INTRODUÇÃO 


O processo de avanço tecnológico nos diversos campos da ciência é perene e 
vem ocorrendo de forma bastante rápida. Uma das consequências disto é o surgimento 
de novos materiais para aplicação em vários ramos da Engenharia. Muitas destas 
aplicações envolvem problemas de transferência de calor, fazendo com que as técnicas 
de caracterização térmica se tornem cada vez mais importantes. O fato de uma boa 
caracterização térmica dos materiais tornar-se um fator crítico de sucesso nos projetos 
de Engenharia justifica o número crescente de pesquisas e desenvolvimento de novas 
técnicas para determinação de propriedades termofísicas. Carollo et al. (2015) atribuem 
a globalização a necessidade de desenvolver novas técnicas para determinar de forma 
rápida, confiável e precisa as propriedades termofísicas dos materiais. 


As propriedades térmicas de um material podem ter valores controversos entre 
fontes de literatura, têm dependência significativa da temperatura e dependem do método 
de determinação (CHIRDON e PATIL, 2011). A determinação de propriedades termofísicas 
não é uma tarefa evidente, sendo os problemas divididos em dois grupos: concepção e 
solução do modelo térmico e montagem experimental (BORGES, 2008). 


Há diversas técnicas para medição das propriedades termofísicas dos materiais. 
Existem aquelas que estimam as propriedades de forma isolada e outras onde a estimativa 
é realizada de forma simultânea. Neste contexto, Carslaw e Jaeger (1959) atribuem às 
técnicas periódicas uma grande importância em medições a baixas temperaturas. Dentre 
as técnicas periódicas, destaca-se a proposta em 1861 pelo físico sueco Anders Jonas 
Angstrom (1814 - 1874), que imprime um fluxo de calor periódico sobre a amostra teste, 
provocando nesta um campo de temperatura também periódico. 


Uma propriedade termofísica extremamente importante é a difusividade térmica. Ela 
reflete a razão entre a energia que determinado meio pode transportar pelo processo de 
difusão e a energia que este pode armazenar, ou seja, a difusividade térmica mostra a 
rapidez com que o calor pode se propagar em determinado material. 


Diante do exposto, esta pesquisa, cuja temática aborda a difusividade térmica de 
determinados materiais, propõe-se a responder o seguinte questionamento: É possível 
estimar a difusividade térmica de determinados materiais metálicos, utilizando um método 
baseado na técnica de Angstrom, e dentro de uma faixa de incerteza aceitável? Para tal 
questão, devem ser admitidas as hipóteses de que, na amostra, o fluxo de calor será 
unidimensional e que não haverá geração de energia. O objetivo deste trabalho é estudar 
teórica e experimentalmente a difusividade térmica de ligas metálicas, limitando-se aos 
aços inox AIS! 304 e AlSI 316. Em termos específicos, esta pesquisa visa: a) Realizar 
a modelagem matemática do problema físico que reflete o objetivo deste trabalho; b) 
Construir um dispositivo experimental baseado na técnica de Angstrom, de modo a viabilizar 
o objetivo deste trabalho; c) Realizar a calibração dos termopares, no intuito de minimizar 
erros nos valores de temperatura a serem medidos; d) Identificar as temperaturas médias, 
amplitudes e fases dos sinais de temperatura captados pelos termopares. 


O presente trabalho foi organizado em 8 capítulos, a contar deste. O capítulo 2 versa 
sobre as equações fundamentais da transferência de calor por condução, a lei de Fourier 
e a equação da difusão de calor, assim como a influência que a propriedade objeto desta 


pesquisa exerce nestas expressões. O capítulo 3 apresenta uma revisão sobre métodos 
de determinação das propriedades termofísicas, com ênfase aos métodos periódicos. O 
capítulo 4 explana a modelagem matemática do problema físico desta pesquisa. O capítulo 
5 discorre sobre o dispositivo experimental construído para assegurar o objetivo deste 
trabalho. O capítulo 6 aborda o procedimento de calibração realizado nos termopares 
soldados nas amostras. No capítulo 7, os resultados encontrados são interpretados, 
analisados, e comparados com valores disponíveis na literatura. No capítulo 8, a conclusão 
do trabalho, explicita-se a resposta à pergunta levantada no capítulo 1, bem como possíveis 
limitações e propostas de futuros trabalhos. 


Introdução 


EQUAÇÕES FUNDAMENTAIS DA CONDUÇÃO DE CALOR E 
PROPRIEDADES TERMOFÍSICAS 


A análise das equações que regem os fenômenos da transferência de calor por 
condução requer o conhecimento de diversas propriedades da matéria, geralmente 
conhecidas como propriedades termofísicas. Estas, por sua vez, encontram-se divididas 
em duas categorias: propriedades de transporte e propriedades termodinâmicas. As 
propriedades de transporte dizem respeito aos coeficientes de taxa de difusão, enquanto 
as propriedades termodinâmicas estão relacionadas ao equilíbrio e mudanças de estado de 
um sistema. Este capítulo trata da influência que as propriedades termofísicas de transporte 
exercem nas equações fundamentais da condução de calor: a lei de Fourier e a equação 
da difusão de calor. 


11 LEI DE FOURIER E A CONDUTIVIDADE TÉRMICA 


De acordo com Ozisik (1993), a lei básica que relaciona o fluxo de calor e o gradiente 
de temperatura, baseada em observações experimentais, ficou conhecida como Lei de 
Fourier após o cientista francês Joseph Fourier a utilizar no seu trabalho, Theorie Analytic 
de la Chaleur, publicada em Paris, no ano de 1822. Incropera e De Witt (2003) enfatizam 
que a Lei de Fourier é a base da transferência de calor por condução e trata-se de uma 
expressão vetorial, indicando que o fluxo de calor é normal a uma superfície isoterma e na 
direção decrescente da temperatura. Aplica-se a todas as substâncias, independentemente 
de seu estado (sólido, líquido ou gasoso). A Lei de Fourier é dada pela Equação 1. 
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onde q” é o fluxo de calor (W/m?), k é a condutividade térmica (W/m.K) e oT/dOn é o 
gradiente de temperatura na direção normal às isotermas (K/m). 


O uso da lei de Fourier só será possível mediante o conhecimento da condutividade 
térmica. Tal propriedade indica a taxa pela qual a energia é transferida pelo processo 
de difusão. A magnitude da condutividade térmica varia conforme o tipo e a composição 
química do material, sua estrutura física e seu estado. Além disso, depende também da 
temperatura do material. A Figura 1 mostra a condutividade térmica como uma função da 
temperatura para várias substâncias. 
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Figura 1 - Condutividade térmica de várias substâncias versus temperatura. 


Fonte: Kakaç e Yener (1993) 


Nota-se, da Figura 1, que o valor da condutividade térmica de um sólido metálico 
possui, em média, uma ordem de grandeza 4 vezes maior do que a dos gases. Os sólidos 
metálicos normalmente são chamados de condutores, por possuírem um valor elevado de 
k, enquanto materiais com baixo k são chamados de isolantes. Ainda conforme a Figura 
1, nota-se que os valores de k para sólidos metálicos, de maneira geral, decrescem com 
o aumento da temperatura. Obviamente, a variação da condutividade térmica com a 
temperatura pode ser negligenciada quando a faixa de temperatura em consideração não 
é tão elevada ou a dependência da condutividade térmica em relação à temperatura não é 
considerável. 


Tratando especificamente de materiais no estado sólido, Incropera e De Witt 
(2003) comentam que, para tal estado da matéria, o transporte de energia térmica é 
composto por dois efeitos: a migração de elétrons livres e através das ondas vibracionais 
da rede cristalina do material. Desta maneira, a condutividade térmica é a soma entre 
o componente eletrônico k, e o componente da rede k. De forma aproximada, k, é 
inversamente proporcional à resistividade elétrica p,. Em metais puros, cujos valores de 
p. São baixos, o valor de k, é muito maior que k.. Já para o caso das ligas metálicas, cujos 
valores de p, são consideravelmente mais elevados, a contribuição de k, para o valor de 
k não pode mais ser desprezada. 
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21 DIFUSIVIDADE TÉRMICA E A EQUAÇÃO DA DIFUSÃO DE CALOR 


A condutividade térmica reflete a quantidade de energia térmica que o meio pode 
transportar através do processo de difusão. Quando o transporte térmico ocorre de forma 
transiente, parte da energia difundida no meio é acumulada internamente na matéria e o 
restante é transferida na direção do gradiente térmico. Nessa condição, outra propriedade 
térmica se torna importante para avaliar o fenômeno, a capacidade calorífica volumétrica. 
Definida pelo produto entre duas propriedades termodinâmicas, massa específica (p) e 
calor específico (C,), seu valor representa a capacidade de um material armazenar energia 
térmica. 


A razão entre a condutividade térmica e a capacidade calorífica volumétrica define 
uma propriedade importante chamada difusividade térmica, a, cuja unidade no SI é ms. 
A difusividade térmica representa a razão entre a capacidade do material conduzir energia 
térmica e sua habilidade de armazená-la. Em outras palavras, ela representa a relação 
entre a energia que o meio pode transportar (representada pela condutividade térmica) 
e a quantidade de energia que ele pode armazenar (representada pela capacidade 
calorífica volumétrica). Materiais com a elevados respondem rapidamente a mudanças nas 
condições térmicas impostas, enquanto materiais com valores reduzidos de a respondem 
mais lentamente, levando um tempo maior para atingir uma nova condição de equilíbrio. 
A difusividade térmica exerce papel fundamental em problemas condutivos no regime 
transiente. 


De maneira geral, o processo de condução de calor é fundamentado na 1º Lei 
da Termodinâmica ou equação da conservação de energia. Considerando, um volume 
de controle infinitesimal num meio homogêneo, sem advecção, onde há um campo de 
temperatura, ou seja, um gradiente de temperatura, que não há geração interna de calor e 
que o meio é isotrópico e com propriedades constantes, a equação da difusão de calor, cuja 
forma geral é dada pela Equação 2!, reduz-se à Equação 3, em coordenadas cartesianas. 


(6) + (65) + (6) + 6 = 9T 

0xN 0x) oyv 0y)' dz2N' dz — Pep'a (2) 
o2T 92T 92T 19T 

e + ef e e cm 3) 
dx? dy? dz? adt ( 


A Equação 3 é uma equação diferencial cuja solução, para condições de contorno 
especificadas, revela a distribuição de temperatura do meio. Sua aparente complexidade 
não deve obscurecer o fato de que ela descreve uma condição física importante, ou seja, a 
conservação da energia. A equação da difusão de calor constitui a base para a solução do 
modelo matemático proposto para a presente pesquisa. 


1 A dedução da equação 2 está disponível em Incropera e De Witt (20083). 
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MÉTODOS DE IDENTIFICAÇÃO DA DIFUSIVIDADE TÉRMICA 


Segundo Carslaw e Jaeger (1959), existe um elevado número de métodos utilizados 
para identificação de propriedades termofísicas. Esses métodos podem ser classificados 
essencialmente em três grupos: permanentes, transientes e de aquecimento periódico, 
estando o presente trabalho situado neste último. Neste capítulo, deu-se um maior destaque 
aos métodos periódicos. Ele subdivide-se em três partes. Na primeira, discutem-se métodos 
considerados clássicos, por conta do vasto uso destes em pesquisas; na segunda, realiza- 
se uma revisão dos métodos periódicos, dentre eles, a técnica de Angstrom, base desta 
pesquisa; na terceira, têm-se outros métodos desenvolvidos em diversos trabalhos, para 
estimativa de propriedades termofísicas. 


11 MÉTODOS CLÁSSICOS 


De acordo com Santos (2005) e Gravena et al. (2010), dois métodos frequentemente 
empregados em laboratórios e centros de pesquisa para determinação de propriedades 
termofísicas são o método do fio quente e o método flash. Nas subseções conseguintes, 
estes métodos dois métodos são brevemente descritos. 


1.1 Método do fio quente 


No tocante ao método do fio quente, pode-se afirmar que se trata de um método 
absoluto, direto e não estacionário (GRAVENA et al., 2010). De acordo com Gravena et 
al. (2010), o procedimento de realização do método do fio quente inicia-se com a inserção 
de uma sonda cilíndrica, no centro axial da amostra que se deseja medir. A sonda tem a 
função de dissipar calor por efeito Joule e medir a temperatura no interior da amostra. 
Teoricamente, o fio quente é considerado como uma fonte de calor ideal, infinitamente 
longo e fino, o qual é circundado infinitamente pelo material cujas propriedades térmicas se 
deseja determinar. Ao passar uma corrente elétrica pelo fio, uma taxa constante de calor 
no tempo e no comprimento da amostra será liberada e se propagará pelo material. Essa 
propagação num meio infinito gera, no material, um campo transiente de temperaturas. 


Conforme Santos et al. (2004), o método do fio quente possui duas limitações: 
materiais condutores elétricos e materiais de alta condutividade térmica. Esta segunda 
limitação também é mencionada como uma limitação do método em Gravena et al. (2010), 
pelo motivo de existir uma alta resistência térmica de contato entre a sonda e a amostra, 
visto que é muito difícil eliminar os interstícios de ar presentes na montagem. 


1.2 Método flash 


O método flash trata-se de um método transiente, cujo princípio de funcionamento 
encontra-se descrito em Santos (2005): “[...] um pulso de energia de curta duração incide 
na face frontal da amostra a ser ensaiada, e a difusividade térmica é calculada a partir 
do registro do histórico de temperatura na face posterior da amostra”. Com o valor da 
difusividade térmica estimado experimentalmente, e conhecendo-se a magnitude do 
calor específico e da massa específica do material, determina-se então a condutividade 
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térmica. O método flash é direto na determinação da difusividade térmica, porém indireto na 
estimativa da condutividade térmica. Maiores detalhes sobre o princípio de funcionamento 
do método flash pode ser verificado em Ryu et al. (2013) e Parker et al. (1961). 


O método flash é largamente utilizado na determinação da difusividade térmica 
de materiais metálicos, cerâmicos, compósitos, metais líquidos, materiais poliméricos 
e determinados alimentos. De acordo com Reif-Acherman (2014), a vantagem mais 
importante associada a este método é a facilidade de preparação da amostra devido ao 
seu tamanho reduzido, além da eliminação do problema da resistência térmica de contato 
e a redução nas perdas de calor. 


2 | MÉTODOS PERIÓDICOS 


De acordo com Carslaw e Jaeger (1959), o conjunto de técnicas denominadas 
técnicas periódicas assumem um papel fundamental na identificação de propriedades 
termofísicas dos materiais quando as medições são realizadas em baixas temperaturas. 
Dentre estas técnicas, destaca-se a que foi desenvolvida, em 1861, pelo físico sueco 
Anders Jonas Angstrom (1814 - 1874). 


Nas subseções a seguir, detalha-se o princípio de funcionamento do método 
desenvolvido por Angstrom, além de outras técnicas periódicas que foram baseadas neste 
método. 


2.1 Método de Angstrom 


Em 1861, Angstrom propôs uma técnica experimental para determinação da 
condutividade térmica do cobre e do ferro. Ao contrário de outros métodos disponíveis àaépoca, 
os quais necessitavam de um regime permanente de temperatura e consequentemente um 
fluxo de calor permanente, no método de Angstrom uma extremidade de uma barra longa 
com seção transversal pequena, porém uniforme, foi sujeita a uma variação periódica na 
temperatura, sendo alternadamente aquecida e resfriada, em intervalos de tempos iguais, 
enquanto a outra extremidade era exposta à temperatura do ambiente. Na Figura 2 verifica- 
se a montagem básica do aparato experimental desenvolvido por Angstrom. 
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Figura 2 - Aparato experimental desenvolvido por Angstrom. 


Fonte - Reif-Acherman (2014) 


Na Figura 2, a válvula indicada por b pode assumir a posição b, ou a posição b,. 
Isto determina quando a amostra, cuja seção transversal identifica-se pôr a, pode ser 
submetida ao vapor quente vindo do gerador de vapor A ou à água resfriada proveniente do 
vaso B. A temperatura flutua periodicamente em cada ponto ao longo da barra. Devido aos 
efeitos da radiação e da convecção na superfície, as ondas térmicas que se deslocam ao 
longo do comprimento da barra são atenuadas, e se propagam com uma diferença de fase. 
Após alguns ciclos, atinge-se o regime periódico, independente das condições iniciais. A 
variação da temperatura ao longo da barra torna-se uma função com a mesma frequência 
da fonte de calor, sendo, então, mensurada em dois pontos ao longo da barra, através de 
um termômetro. 


2.2 Métodos periódicos fundamentados na técnica de Angstrom 


Em Haji-Sheikh et al. (1998) reporta-se uma técnica baseada no método de 
Angstrom para determinação simultânea da condutividade térmica e da difusividade 
térmica. Inicialmente realizou-se uma análise de sensibilidade para determinar-se a melhor 
faixa de frequência para este tipo de investigação experimental. Os materiais utilizados 
foram o Delrim [um tipo de polímero] e o aço inox AlISI 304. Os resultados obtidos para 
a condutividade térmica do aço inox AlISI 304 e para o Delrim, quando comparados aos 
valores disponíveis na literatura, tiveram erros de 9% e 20%, respectivamente. 

Bezerra Filho et al. (1999) apresentaram uma técnica que utiliza O sinal periódico 
para determinar a difusividade térmica. Um dispositivo experimental foi construído, a fim 
de realizar medidas precisas da difusividade térmica utilizando um sinal térmico periódico 
de frequência elevada, entre 0,01 Hz e 0,25 Hz. A amostra utilizada foi feita de cobre, com 
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42 mm de diâmetro e 100 mm de altura. A face da amostra em contato com o aquecedor 
ficou submetida a uma fonte de calor periódica, enquanto isso, a face inferior foi mantida a 
uma temperatura aproximadamente constante, mediante interação com um fluido térmico. 
A evolução da temperatura na amostra foi acompanhada por 7 termopares nela soldados. 
Estes termopares são do tipo K com diâmetro de 80 um e foram soldados através de 
descarga capacitiva. Dessa forma, a difusividade térmica foi obtida, ora utilizando-se a 
razão das amplitudes dos sinais térmicos obtidos pelos termopares, ora utilizando-se a 
defasagem com relação a um ponto de referência. Devido ao vácuo gerado pelo dispositivo, 
o número de Biot é muito menor que 1 e, dessa forma, a temperatura em cada seção 
transversal da peça é uniforme. O modelo matemático pôde então ser formalizado como 
unidimensional. Os resultados obtidos por Bezerra Filho et al. (1999) quando comparados 
com resultados obtidos através do método flash para o mesmo material apresentaram um 
desvio de 15%. 


Lahoucine e Khellaf (2004) utilizaram um método periódico para determinação 
simultânea e precisa da difusividade e da condutividade térmicas, considerando nos cálculos 
os efeitos da presença dos termopares. A difusividade térmica foi obtida similarmente ao 
trabalho de Haji-Sheikh et al. (1998), com o uso da defasagem do sinal térmico. Por outro 
lado, na determinação da condutividade térmica foram considerados fatores de correção 
devido aos efeitos gerados pelos termopares. Essa condição foi negligenciada por Haji- 
Sheikh et al. (1998) e que, segundo Lahoucine e Khellaf (2004), pode ter sido a causa dos 
elevados erros experimentais encontrados naquele trabalho. A presença dos termopares 
provoca distúrbios localizados em torno do ponto de medição, cuja importância destes 
distúrbios vai depender tanto da natureza do contato térmico quanto das propriedades 
termofísicas do sensor de temperatura e do material sólido onde os termopares estão 
fixados. A solução do modelo matemático aplicado a este método foi obtida através do 
teorema de Duhamel. 


Santos et al. (2010) utilizaram uma variação simples e de baixo custo do método 
de Angstrom em procedimentos experimentais para determinação da difusividade térmica 
de determinados polímeros. Os resultados alcançados através destes procedimentos 
atingiram uma boa precisão, entre 0,1 e 1%, quando comparados aos valores identificados 
através dos métodos do fio-quente e flash. 
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O MODELO MATEMÁTICO 


Na Figura 3 pode ser visualizado o sistema físico que representa o problema 
apresentado neste trabalho. 


T(L,t) = A(L).sen(wt + E) 


Amostra 





Figura 3 - Sistema físico do problema. 


Fonte - Autoria própria (2017) 


Para chegar-se ao modelo matemático que representa o sistema físico da Figura 3 
admite-se que os experimentos seriam ora realizados no vácuo, reduzindo desta maneira o 
coeficiente de convecção e consequentemente aumentado a resistência térmica convectiva 
na direção radial, e ora com o uso de isolante térmico em contato com a área lateral da 
amostra, aumentando assim a resistência térmica radial, garantindo que o fluxo de calor 
ocorresse na direção axial. Desta feita, tem-se que o número de Biot [cuja definição pode 
ser encontrada em Incropera e De Witt (2008, p.172)] assumiria valores menores que 0,1 
para as condições experimentais em questão, implicando, portanto, que a temperatura 
em cada seção da amostra seria uniforme durante o transporte de energia, condição para 
que a transferência de calor por condução, ao longo da amostra, ocorresse de maneira 
unidimensional. Supõe-se ainda que a variação da difusividade térmica com a temperatura 
é desprezível e que o meio é isotrópico com propriedades constantes. Isto posto, 
considerando uma amostra de comprimento L (m), a Equação 3 reduz-se à Equação (4.1). 


o2T 107 
dx? adt (3) 


onde T(x,t) é o campo térmico (ºC), x é a componente espacial (m), a é a difusividade 
térmica (m?/s) e t é o tempo (s), com O <x<Let> 0. A Equação (4.1) é uma equação 
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diferencial parcial de segunda ordem e, para resolvê-la, faz-se necessário uma condição 
inicial e duas condições de contorno. A condição inicial que reflete o problema físico é 
dada pela Equação (4.2.a). As duas condições de contorno foram obtidas nas extremidades 
da amostra. Sobre a parte superior, foi imposto um fluxo de calor periódico com uma 
dada frequência. A extremidade inferior foi mantida a uma temperatura constante, por 
ter ficado em contato com um fluido cuja temperatura é controlada através de um banho 
termorregulável. As expressões que representam as condições nas faces superior e inferior 
da amostra são dadas respectivamente pelas Equações (4.2.b) e (4.2.0). 


Metj=0 QUasval (4.2.3) 
T(L,t) = A(x,) sin(wt + €) (4.2.b) 
T(O,t) =0 (4.2.C) 


onde A(x,) = 1 é a amplitude do modelo na posição x = L, e w é a frequência do sinal 

térmico (rad/s). A Equação 5 mostra o campo de temperatura solução deste modelo, na 
forma adimensional, obtida por Carslaw e Jaeger (1959, p.105). 

n(=1)" [an2n? sen(e) — wL? cos(£)] ad (=) (Ee) 


T(x, tb) = t W++2 >. 
(x,t) sen(wt+e+V) ++27a ninint + ilá 


n=1 


(5) 


onde T é o campo de temperatura adimensional, A e ly são respectivamente a 
amplitude e a fase, do modelo, do sinal térmico captado em determinado ponto ao longo 
da amostra, e a é a difusividade térmica (m?/s). As expressões que definem A e IP são 
dadas, respectivamente, pelas Equações 5.1 e 5.1. A dedução destas expressões pode ser 
verificada no Apêndice A. 


senh[xp (1 + D]| Eesti — cos(28%)] 2 








“IsenhLB +D]|  I|cosh(28L) — cos(2BL) Ee 
o senh[xf (1 + i)] 
P= arg ET] (5.2) 
Onde 
— [o 
“Za 


Na Equação 5, a perturbação transiente ligada à condição inicial tende para zero, 
conforme mostrado por Carslaw e Jaeger (1959), e o campo de temperatura passa a ser 
permanente periódico. Pode-se obter, de forma alternativa, uma solução do campo de 
temperatura permanente periódico quando este regime for alcançado. Para isto, considera- 
se que as variações transitórias do campo de temperatura cessam quando o tempo cresce 
(o termo transitório da Equação 5 desaparece quando t — ºº) de tal maneira que em tempos 
longos atinge-se a condição de regime periódico permanente. Esta solução é dada pela 
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Equação 5. 


T(x,t) = Asen(wt + e +'P) (5) 


onde À, w, €, e IP já foram definidos anteriormente. 
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DISPOSITIVO EXPERIMENTAL 


Este capítulo trata dos sistemas e seus componentes, que constituem o dispositivo 
experimental usado na pesquisa além, dos procedimentos de operação deste dispositivo. 


O dispositivo experimental ilustrado na Figura 4, concebido para identificação da 
difusividade térmica utilizando um campo de temperatura periódico, é constituído pelos 
sistemas de aquisição de dados, de aquecimento, de resfriamento, de vácuo, do cabeçote 
de ensaio e da amostra. 


Hardware 

Gaiola de Faraday 
Cabos de comando e controle 
Aquisitor de dados 
Termopares tipo K 

Fonte de potência 
Cabeçote de ensaio 
Saída de fluido térmico 
Retorno de fluido térmico 
Banho termorregulador 
Linhas de vácuo 


Vacuômetro 


A 
B 
c 
D 
E 

F 

G 
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M 


Bomba de vácuo 





Figura 4 - Dispositivo Experimental. 


Fonte - Autoria própria (2017) 


11 SISTEMA DE AQUISIÇÃO DE DADOS 


O papel fundamental do sistema de aquisição de dados é a coleta e o armazenamento 
dos dados dos sinais de temperatura gerados no decorrer dos experimentos. Este sistema 
é constituído por um aquisitor de dados 34970 da Agilent, o qual possui um componente 
interno que dispensa a utilização da chamada junção de referência ou junção fria no uso de 
termopares, e por um microcomputador, no qual há uma rotina, em Visual Basic, que realiza 
a aquisição dos dados e estabelece o número total de ciclos do experimento 
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21 SISTEMA DE AQUECIMENTO 


O sistema de aquecimento compreende uma resistência elétrica de 0,26 O e 
uma fonte de potência E3633A da Agilent, comandada pela mesma rotina que controla a 
aquisição de dados. O papel da fonte de potência é estabelecer um fluxo de calor periódico 
para a amostra. Entre outras palavras, o objetivo da fonte é produzir um sinal de tensão 
retificado na resistência elétrica, conforme Equação 7. 


V = Yin [sen (5º) (7) 


onde V é a tensão gerada pela fonte (volts), V. é a tensão média (volts), w é a 
frequência térmica (rad/s) comandada pelo sistema que controla o aquecimento etéo 
tempo (s). Desta maneira, a resistência elétrica gera uma potência dada pela Equação 8. 


P = B,a[1 — cos(wt)] (8) 


onde P é a potência gerada (W) e P. é a potência média (W). 


31 SISTEMA DE RESFRIAMENTO 


O sistema de resfriamento é composto pelo fluido térmico 47V10 da Siliplus e por um 
banho termorregulador 12101-56 da Cole-Parmer. Este possui internamente um reservatório 
onde o fluido térmico é armazenado e onde há uma bomba centrífuga, responsável pela 
circulação do fluido, entre o banho termorregulador e o cabeçote de ensaio. A função deste 
sistema é garantir que, durante a realização do experimento, a temperatura da face inferior 
da amostra esteja à 0ºC, condição estabelecida no modelo matemático. 


41 SISTEMA DE VÁCUO 


O sistema de vácuo é composto por uma bomba de vácuo Speedvac 2 da 
Edwards, com capacidade de gerar uma pressão negativa de 710 mmHg no sistema, e um 
vacuômetro com coluna de mercúrio. A função deste sistema é gerar vácuo no interior da 
câmara da amostra numa intensidade suficiente para manter o número de Biot neste local 
menor que 0,1, condição necessária para assegurar a hipótese de transferência de calor 
por condução unidimensional, adotada na solução do modelo matemático. Este sistema 
deve ser desativado em experimentos que utilizarem isolante térmico ao invés do vácuo. 


5 | AMOSTRAS 


Foram estudadas duas amostras de formato cilíndrico, constituídas de três materiais 
diferentes, respectivamente, aço inox AlSI 304 (amostra 1) e aço inox AlSI 316 (amostra 2). 
Os processos de fabricação das amostras e o processo de instalação dos termopares são 
descritos nas subseções conseguintes. 
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5.1 Fabricação das amostras 


As amostras de aço inox AISI 304 e aço inox AlSI 316, respectivamente, foram 
confeccionadas pelo processo de torneamento, a partir de tarugos comerciais trefilados 
com 12,7 mm de diâmetro e 170 mm de altura, de tal maneira que a dimensão definitiva de 
cada amostra, após o torneamento, foi de 12,7 mm de diâmetro e 150 mm de altura. Em 
seguida, as faces superiores e inferiores foram lixadas. A face superior de cada amostra 
também foi polida, com o objetivo de reduzir a RTC (Resistência Térmica de Contato) entre 
esta face e a resistência elétrica. 


5.2 Instalação dos termopares 


Em cada amostra deste trabalho foram instalados 14 termopares do tipo K (Cromel- 
Alumel) com 100 um de diâmetro e 1 m de comprimento. A opção por este tipo de termopar 
deve-se ao fato deste apresentar um comportamento linear em sua curva característica 
tensão versus temperatura numa ampla faixa, quando comparado aos outros tipos de 
termopares. 


A fixação dos termopares nas amostras se deu através do processo de soldagem por 
descarga capacitiva. Esta técnica tem a vantagem de assegurar um contato perfeito entre 
Os termopares e os pontos onde deseja-se medir as temperaturas e de reduzir o tempo de 
resposta que poderia provir de uma resistência térmica de contato entre os termopares e o 
corpo-de-prova (BEZERRA FILHO et al., 1999). A Figura 5 ilustra o processo de soldagem 
dos termopares. 


Amostra 


Detalhe da soldagem do termopar tipo K à amostra 


Eletrodo da solda por descarga capacitiva 
Fio CROMEL (+) 

Fio ALUMEL (-) 

Cabo de termopar tipo K 





Figura 5 Soldagem de Termopar. 


Fonte - Autoria própria (2017) 
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5.2.1 | Identificação e posição de soldagem dos termopares na amostra 


Buscando acelerar a montagem dos termopares no aquisitor de dados, além de 
facilitar a apresentação dos resultados das razões de amplitude e defasagem entre os 
sinais de temperatura, é adotada uma nomenclatura sequencial para identificação dos 
termopares, conforme visualizado na Figura 6. 





Figura 6 - Identificação dos termopares. 


Fonte - Autoria própria (2017) 


Como pode ser visto na Figura 6, foram soldados 14 termopares em cada amostra. 
Para identificação da razão de amplitudes e da defasagem e, consequentemente, da 
difusividade térmica, foram utilizados os termopares “xr” e os termopares “x1” até “x5”, 
sendo o termopar “x6” utilizado apenas para observar o comportamento da atenuação 
da amplitude das oscilações periódicas de temperatura. Os termopares marcados com 
(*) são os termopares de reserva, e seriam utilizados caso o respectivo termopar titular 
apresentasse algum defeito durante o experimento. As distâncias dos locais de fixação 
dos termopares, ditos titulares, para as faces da amostra é a mesma para os respectivos 
termopares de reserva. 


Vale salientar que, para evitar possíveis efeitos bidimensionais de transferência de 
calor decorrente da proximidade da resistência elétrica, o termopar de referência (xr) e seu 
substituto (xr*) foram instalados a uma certa distância daquele componente. Com isso, o 
local onde o termopar de referência está instalado foi considerado como sendo o final da 
amostra. O termopar de referência nas amostras de aço inox AlSI 304 e de aço inox AlSI 
316 foi instalado a uma distância de 20 mm da resistência elétrica. 


Com isso, nos modelos matemáticos utilizados para identificar A, VW e a, o 
comprimento das duas primeiras amostras foi considerado como sendo 130 mm. A Tabela 
1 apresenta a posição de cada termopar em relação ao termopar de referência, para as 
amostras em estudo neste trabalho. 
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Termopar Amostra 1 Amostra 2 


xr 130 130 
x1 120 110 
x2 115 105 
x3 110 100 
x4 105 95 
x5 100 90 
x6 5 5 


Tabela 1 - Posição (mm) dos termopares em relação ao termopar de referência “xr”. 


Fonte - Autoria própria (2017) 


61 CABEÇOTE DE ENSAIO 


O cabeçote de ensaio é um dispositivo em formato cilíndrico. É o local onde a 
amostra fica durante o experimento. Pode ser visto em detalhes na Figura 7. 


Reservatório de fluido 
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Gaiola de Faraday 





Figura 7 - Cabeçote de ensaio. 


Fonte - Autoria própria (2017) 
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Conforme a Figura 7, os componentes cujos balões de identificação estão na cor azul 
são aqueles que integram o cabeçote de ensaio, dentre os quais destaca-se o componente 
B, denominado câmara da amostra. 


A câmara da amostra é constituída por um tubo de parede fina de aço inox com 123 
mm de altura e cujo valor do diâmetro, 100 mm, foi dimensionado a partir do conceito de 
raio crítico de isolamento. Isto se deve ao fato de existir a possibilidade do experimento 
ser realizado ou com vácuo, ou com isolante térmico (neste trabalho utilizou-se isopor) no 
interior da câmara. Desta maneira, o diâmetro dimensionado foi aquele para o qual a perda 
de calor radial pelo isolamento térmico apresentou um valor desprezível. 


Ainda compõem a câmara da amostra duas tampas para fechamento, feitas em 
PVC, que possuem um furo centralizado de 12,7 mm, mesmo diâmetro da amostra. Na 
tampa superior há mais dois furos: um com 3 mm de diâmetro, por onde serão passados 
os termopares e outro com 10 mm de diâmetro, para instalação de uma conexão para 
realização de vácuo no interior da câmara. Na tampa inferior, em contato com a face inferior 
da amostra e com o fluido de refrigeração, há um disco de cobre de 2 mm de espessura, 
instalado para evitar a entrada deste fluido na câmara da amostra por capilaridade. Admitiu- 
se que este disco, por conta da alta condutividade térmica do cobre, é isotérmico em 
relação ao fluido de refrigeração. A câmara com a amostra em seu interior, é montada sobre 
o reservatório de fluido, local por onde circula o fluido de refrigeração, a 0ºC. 


Vale destacar, que o cabeçote de ensaio, o aquisitor de dados e a fonte de potência, 
durante os experimentos, ficaram no interior de uma gaiola metálica, baseada no princípio 
da Gaiola de Faraday. Esta gaiola foi utilizada para minimizar os efeitos das perturbações 
causadas pelas interferências eletromagnéticas decorrentes de correntes parasitas, 
conhecidas como Correntes de Foucault. 
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CALIBRAÇÃO DOS TERMOPARES 


No intuito de minimizar os erros dos valores de temperatura a serem medidos, 
realizou-se procedimentos de calibração dos termopares instalados nas amostras. Estes 
procedimentos consistem no levantamento de curvas de calibração, as quais foram utilizadas 
para corrigir os valores de temperatura obtidos [nos experimentos para identificação 
de a] em valores reais de temperatura. Neste capítulo, detalha-se como os termopares 
foram calibrados. Além disto, apresentam-se as curvas de calibração levantadas para os 
termopares das duas amostras utilizadas. 


11 PROCEDIMENTOS DE CALIBRAÇÃO 


A curva da tensão termoelétrica em função da temperatura para o tipo de termopar 
utilizado neste trabalho, tipo K, apresenta um comportamento linear para uma ampla faixa 
de temperatura (-200 até 1250ºC). Diante disto, foi necessário o levantamento de dois 
pontos, no plano coordenado, para obtenção da curva de calibração de cada termopar, em 
cada amostra. Estes pontos foram as temperaturas de fusão do gelo e de ebulição da água 
para a cidade de Campina Grande, Paraíba, Brasil. Vale salientar que a água utilizada em 
ambos os casos foi do tipo destilada. Para cada uma das três amostras foram realizados, 
então, dois procedimentos de calibração. 


A temperatura de fusão do gelo, para a referida cidade é de aproximadamente 0,0ºC, 
valor obtido a partir da equação de Clapeyron para o equilíbrio de fase sólido-líquido. Por 
sua vez, a temperatura de ebulição da água, por ser dependente da pressão atmosférica, 
foi obtida da seguinte maneira: por meio de um barômetro, aferiu-se a pressão atmosférica 
no instante da realização do procedimento de calibração. Esta pressão foi inserida no 
software CATT (Computer Aided Thermodynamic Tables), onde obteve-se a temperatura 
de vapor saturado da água em função de tais condições de pressão. Para o caso das 2 
amostras em análise, essa temperatura foi de aproximadamente 98,0ºC. 


Vale salientar que, em ambos os procedimentos de calibração, no ponto de fusão do 
gelo e no ponto de ebulição da água, as amostras ficaram na posição horizontal, evitando 
que possíveis efeitos de convecção natural interferissem nas medidas dos termopares 


Ao final dos procedimentos de calibração, foram plotadas as curvas de calibração, 
de cada termopar para cada amostra. As equações destas curvas são listadas na Tabela 2. 


Calibração dos termopares EEN 


Amostra 
Aço inox AISI 304 
Aço inox AISI 304 
Aço inox AISI 304 
Aço inox AlSI 316 
Aço inox AISI 316 
Aço inox AlSI 316 


Termopar 
xr e 
Xl, Xe sent 
X3, x4, xD J 
xr, x4 de 
XI; X2,%9 E aBii 
x5 T, 


REAL T 


REAL 


REAL 


Curva de calibração 

-0,8091 + 1,0114 x T 
= -0,9102 + 1,0114 x T 
-0,9093 + 1,0103 x T 
=-0,9121 + 1,0134 x T 
=-0,9112 + 1,0124 x T 
= -0,8099 + 1,0124 x T 


MEDIDA 


MEDIDA 


MEDIDA 


MEDIDA 


MEDIDA 


MEDIDA 


Tabela 2 - Curvas de calibração dos termopares para cada amostra. 


Fonte - Autoria própria (2017) 


Calibração dos termopares Bo. 


RESULTADOS E DISCUSSÕES 


Foram realizados 9 experimentos, visando identificar os valores de a (difusividade 
térmica) das amostras estudadas, a partir dos valores de A (amplitude do modelo ou razão 
de amplitudes) ou de !Y (fase do modelo ou defasagem) dos sinais térmicos captados pelos 
termopares instalados nas amostras. 


Na amostra de aço inox AlSI 304 realizou-se seis experimentos, sendo que em três 
deles utilizou-se vácuo na câmara da amostra e nos demais utilizou-se isolante térmico. Na 
amostra de aço inox AlISI 316 realizou-se três experimentos, todos com isolante térmico na 
câmara da amostra. 


Os valores de w adotados para os experimentos nas amostras de aço inox AlSI 
304 e de aço inox AISI 316 foram 4,5 x 103, 8,5 x 10% e 15,3 x 10% rad/s, sendo estes 
equivalentes, em termos de período, a 1396,3, 739,2 e 410,7 s, respectivamente. Na Figura 
8 pode-se visualizar o resultado do perfil de temperatura, decorrente do fluxo de calor 
periódico, para o experimento realizado na amostra de aço inox AISI 304, com w = 8,5 x 
10% rad/s. 


Comprova-se desta maneira, o que foi mostrado por Carslaw e Jaeger (1959), para 
o campo de temperatura resultante numa amostra submetida à um fluxo de calor periódico: 
a medida que o tempo aumenta, a perturbação transiente é dissipada, e o campo térmico 
torna-se um campo com oscilações periódicas permanentes. Constitui-se, dessa maneira, 
o “regime periódico estabelecido”. 


Transiente ligado à condição inicial Regime periódico permanente 





Termopar "xr" 
Termopar "x1" 
Termopar "x2" 
Termopar "x3" 
Termopar "x4" 
Termopar "x5" 
Termopar "x6" 


Variação das Oscilações 

oscilações periódicas 

periódicas permanentes 

devido à não influenciadas 

condição inicial pela condição 
inicial 


pm 
O) 
(o) 
Sm 
(U 
Am 
ED) 
ra 
(U 
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D 
O 
= 
LD) 
t- 
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Figura 8 - Perfil de temperatura na amostra de aço inox AlSI 304 submetida a um fluxo de calor 
periódico. 


Fonte - Autoria própria (2017) 
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A amplitude térmica registrada pelo termopar “x6” — termopar que fica mais próximo 
da face inferior da amostra, e consequentemente do fluido refrigerante — é praticamente 
nula no regime periódico estabelecido. Isto mostra que, a onda térmica gerada pelo fluxo de 
calor periódico dissipa-se ao longo do comprimento da amostra, à medida que se aproxima 
da face inferior desta. Os sinais de temperatura registrados pelo termopar “x6” não foram 
utilizados para determinação de A e UU. 


Os valores de A e IV foram identificados a partir dos sinais de temperatura obtidos no 
regime periódico estabelecido, ou seja, obtém-se estes parâmetros a partir do instante no 
qual a perturbação transiente ligada à condição inicial encontra-se dissipada. Tal instante 
foi denominado como sendo o instante inicial (t = O s) do regime periódico estabelecido. 
Da Figura 9 até a Figura 17 podem ser vistas as curvas ajustadas (através do algoritmo de 
Levenberg-Marquardt) dos perfis de temperatura no regime periódico estabelecido (RPE), 
para todos os experimentos realizados. 
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Figura 9 - Regime periódico permanente: aço inox AISI 304; w = 4,5 x 10% rad/s (vácuo). 


Fonte: Autoria própria (2017) 
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Figura 10 - Regime periódico permanente: aço inox AISI 304; w = 8,5 x 10% rad/s (vácuo). 


Fonte - Autoria própria (2017) 
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Figura 11 - Regime periódico permanente: aço inox AISI 304; w = 15,3 x 10% rad/s (vácuo). 


Fonte - Autoria própria (2017) 
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Figura 12 - Regime periódico permanente: aço inox AISI 304; w = 4,5 x 10% rad/s. 


Fonte - Autoria própria (2017) 
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Figura 13 - Regime periódico permanente: aço inox AISI 304; w = 8,5 x 10% rad/s. 


Fonte - Autoria própria (2017) 
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Figura 14 - Regime periódico permanente: aço inox AISI 304; w = 15,3 x 10% rad/s. 


Fonte - Autoria própria (2017) 
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Figura 15 - Regime periódico permanente: aço inox AISI 316; w = 4,5 x 10% rad/s. 


Fonte - Autoria própria (2017) 
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Figura 16 - Regime periódico permanente: aço inox AISI 316; w = 8,5 x 10% rad/s. 


Fonte - Autoria própria (2017) 
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Figura 17 - Regime periódico permanente: aço inox AISI 316; w = 15,3 x 10% rad/s. 


Fonte - Autoria própria (2017) 
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Das Figuras 9 - 17, identificam-se A(x), amplitude térmica, W (x), fase, e T (x), 
temperatura média, através dos perfis de temperatura. Estipula-se T(x.t), A(x) e W(x) a 
partir das Equações 9. 


T(x,t) = Tn(x) + Aj(x) senfwt + "'P;(x)] (9.1) 
A(x) = Aj(x)/A;(x,) (9.2) 
W(x) = "Pi(x) — £ (9.3) 


As Figuras 18- 26 apresentam a temperatura média, a razão de amplitudes fases da 
temperatura, em função do posicionamento do termopar. 
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Figura 18 - Temperatura média vs posição do termopar para o aço inox AlSI 304 (vácuo). 


Fonte - Autoria própria (2017) 
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Figura 19 - Temperatura média vs posição do termopar para o aço inox AlSI 304. 


Fonte - Autoria própria (2017) 
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Figura 20 - Temperatura média vs posição do termopar para o aço inox AlSI 316. 


Fonte - Autoria própria (2017) 
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Figura 21 - Razão de amplitudes vs frequência térmica para o aço inox AlSI 304 (vácuo). 


Fonte - Autoria própria (2017) 
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Figura 22 - Razão de amplitudes vs frequência térmica para o aço inox AISI 304. 


Fonte - Autoria própria (2017) 
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Figura 23 - Razão de amplitudes vs frequência térmica para o aço inox AlSI 316. 


Fonte - Autoria própria (2017) 
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Figura 24 - Defasagem vs frequência térmica para o aço inox AISI 304 (vácuo). 


Fonte - Autoria própria (2017) 
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Figura 25 - Defasagem vs frequência térmica para o aço inox AISI 304. 


Fonte - Autoria própria (2017) 
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Figura 26 - Defasagem vs frequência térmica para o aço inox AlISI 316. 


Fonte - Autoria própria (2017) 


No tocante aos gráficos de T| versus posição, verifica-se um comportamento 
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linear crescente de T, conforme o termopar se distancia da origem, com excelente ajuste 
de curva. Isto mostra que os valores considerados realmente encontram-se no regime 
periódico permanente. 


Quanto aos gráficos de A (amplitude do modelo) em função de w, verifica-se que os 
valores de A diminuem de forma exponencial com o aumento de w. Por sua vez, os gráficos 
de !Y (fase do modelo) em função de w, indicam que os valores de IY aumentam de maneira 
exponencial conforme cresce o valor de w. 


Para cada experimento, foram identificados 5 valores de a através da Equação 5.1 e 
5 valores através da Equação 5.2. A partir destes valores, conforme o Guide to Expression 
of Uncertainty in Measurement, do BIPM (2008), foi obtido o valor médio de a e a incerteza- 
padrão u(a). As Figuras 27 - 29 mostram os valores médios de a identificados para cada 
material em função de w, assim como as incertezas-padrão, na parte superior de cada 
amostra. 
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Figura 27 - Difusividade térmica do aço inox AISI 304 (vácuo). 


Fonte - Autoria própria (2017) 
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Figura 28 - Difusividade térmica do aço inox AlSI 304. 
Fonte - Autoria própria (2017) 
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Figura 29 - Difusividade térmica do aço inox AlSI 316. 
Fonte - Autoria própria (2017) 


A partir das figuras 27 - 29, constata-se que o isolante térmico é mais eficiente que 
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o vácuo nas condições em que os experimentos foram realizados. Vale salientar que os 
experimentos no aço inox AlSI 316 foram executados apenas com o uso de isolante térmico 


na câmara da amostra. 


Dentre os valores médios de a identificados nos aços inox AISI 304 e AISI 316, 
através da amplitude ou da fase, selecionou-se aqueles que mais se aproximaram dos 
valores de Carollo et al. (2012). Desta maneira, foram realizados estudos comparativos, 


que podem ser visualizados na Tabela 3. 


Material 


Aço inox 304 
Aço inox 316 


[a + u(a)] a x 10º (m?/s) É á 

x108 (m2/s) Carollo etal. (2012) Sa o Ty (76) 
Ste 0,06 E A -11,41 -1,33 
3,36 + 0,02 3,46 10,12 -2,89 


Tabela 3 - Comparação de valores de difusividade térmica dos aços inox AISI 304 e AISI 316. 


Fonte - Autoria própria (2017) 


De acordo com a Tabela 3, pode-se afirmar que os valores médios de a do presente 
trabalho, para os aços inox AlSI 304 e AISI 316 apresentam desvios cerca de 10% quando 
comparados aos obtidos por Carollo et al. (2012). 

A difusividade térmica trata-se de uma propriedade termofísica que depende da 
temperatura do material. Neste trabalho, pode-se desprezar a variação de a para pequenas 
variações de temperatura, haja visto que os valores de máxima variação de temperatura 
ocorrida em cada amostra são pequenos. 
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CONCLUSÕES 


Este trabalho destinou-se à identificação da difusividade térmica do aço inox AISI 
304 e do aço inox AISI 316, utilizando um campo de temperatura periódico. Para viabilizar 
tal estudo, utilizou-se um método baseado na técnica de Angstrom, na qual faz-se o uso 
de um fluxo de calor periódico na amostra, ocasionando nesta um campo de temperatura 
periódico. Um dispositivo experimental foi desenvolvido para suportar este método. Os 
parâmetros A (razão de amplitudes) e lW (defasagem) entre os perfis de temperatura 
medidos pelos termopares, quando determinados, foram inseridos em seus respectivos 
modelos matemáticos propiciando desta forma a identificação da difusividade térmica dos 
materiais em questão. 


Os resultados obtidos, quando confrontados à literatura, podem ser considerados 
de boa precisão, tendo em vista que as diferenças percentuais ficaram cerca de 10%. Este 
patamar de precisão corrobora com a hipótese de fluxo de calor unidimensional, adotada 
para resolução do modelo matemático deste trabalho. 


Como sugestões para trabalhos futuros, lista-se: 


a. Determinar a condutividade térmica utilizando o mesmo dispositivo experimen- 
tal e/ou através do DSC (neste caso identifica-se o C, e, através da definição da 
a, identifica-se o valor de k). 


b. Estimar a difusividade térmica [utilizando um campo de temperatura periódico] 
de ligas de memória de forma produzidas no LaMMEA cujas temperaturas de 
transformação de fase estejam dentro da faixa de trabalho do dispositivo expe- 


rimental. 


c. Identificar a difusividade térmica [utilizando um campo de temperatura periódi- 


co] de polímeros. 


d. Submeter as amostras deste trabalho aos outros métodos de determinação de 
propriedade termofísicas, como o método flash, e comparar os resultados aos 


obtidos na presente pesquisa. 
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APÊNDICE A - SOLUÇÃO DO MODELO MATEMÁTICO 


O modelo matemático do problema desta pesquisa, a condição inicial e as condições 
de contorno são, respectivamente: 


9ºT OT 
“9x2 dt 
T(x,0) = 0 
T(,)=0 


T(L,t) = sen(wt + €) 


As oscilações transitórias ligadas à condição inicial cessam quando o tempo cresce, 
como mostrado no gráfico da Figura 8 e também mostrado por Carslaw e Jaeger (1959). 
Para tempos longos, só as oscilações permanentes (“regime periódico permanente”) 
continuam. Quando o regime periódico permanente é atingido, o modelo passa a ser: 


927 OT 
“9x2 dt 
T(,0)=0 


T(L,t) = sen(wt +) 


Para resolução deste modelo, segue-se o método apresentado por Myers (1998), 
que consiste em escrever um modelo defasado de 90º. Chamando a variável dependente 
de w no modelo defasado de 90º, tem-se: 


92w —-0w 
“ox dt 
w(O,t) = 0 


w(L,t) = cos(wt + €) 
A partir destes dois modelos, constrói-se o modelo: 
O(x,t) = w(x,t) +iT(x,t) 
Desta forma, o novo modelo é dado por: 
920 00 
md ue (A.1) 
Ox? Ot 
o(0,t) =w(O,t) +IiT(O,t) = 0 
O(L,t) = w(L,t) +iT(L,t) = cos(wt + £) + isen(wt + €) = eitut+e) 
De acordo com Carslaw e Jaeger (1959) e Myers (1998), a solução do modelo é 
dada por: 
9(x,t) = ue Lt = u(y)[cos(wt + €) + isen(wt + €)] (A.2) 


A parte real da Equação (A.2) é a solução do problema w(x,t) e a parte imaginária 
é a solução do problema T(x,t). Substituindo a Equação (A.1) na Equação (A.2), tem-se: 


Apêndice A - Solução do modelo matemático EE 


d0 du 


pad — eilwt+e) 
dx dx 
2 2 
0“0 = õ“u i(wt+e) 
dx? dx? 
do — uimei(wt+e) 
ot 
Com isso, 
du | du q : 
—— > —. — |— — 
x o uiw o ) E u 


Retornando para a Equação (A.2), tem-se: 


(0,0) =u(OMetwv+o =0>u(0)=0 
9(L, t) x u(L)eitwt+e) = eilwt+e) EN u(L) -= 1 


Assim, determina-se u(x) a partir da solução de: 


du wo A (A.3) 
dx? "qr 

u(0) = 0 

= 


A equação característica da Equação (A.3) é: 


2 Ro Ro 
sé-i-=0>s=+|i— 
x 


= [12=0+08 
s=-fiÉ=0+08 


onde: 
ju 
oa 
Assim: 


u(x) = c,e3*% + c,e5* = cyeU+iBx + ce +DBx 


Como u(0) = 0, tem-se: 


u(0) =c+c)=0>c,=-—c 


De u(1) = 1, resulta: 


u(L) = cjeU+DÊL 4 ce U+DBL = q 
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Como c =-c,, resulta: 
cjeG+DBL — ce U+DBL = q (A.4) 
cifeL+DBL — ce-M+DBL] = 1 

Por definição, tem-se: 


(e* —- e*) 
2 
Pode-se escrever a Equação (A.4) como 


senh(x) = 


c;2senh[(1 + DBL]=1 
E, desta maneira: 
1 
“2 2senh[(1 + DPL] 
a 
2 senh[(1 + DL] 


Logo, u(x) é dada por: 


(Q=—(5C0,= 


u(x) = cyel+Bx + ce U+DBx 


e+DBx o(1+DBx 


“0 = eenh[C + DPL] Zsenh(L+ DAL] 
- 2senh[(1 + Dx] 
“CO = Seenh[Cd + DALI 
senh[(1 + )Bx] 
uix) = 


“ senh[(1+ DPL] 


Desta maneira: 


9(x,t) = u(x)eitwt+e) 
= senhiC1 + 262] i(wt+e) 
Blat) — [Eni + | 


O número complexo sm tem uma amplitude (A) e uma fase (lb). Logo, tem- 
senh[(1+i)fL] 





se que: 
a + bi 
u(x) cx di A(x)e 
onde: 
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o a2 + b? 
— qce+d? 
VW=0,-0, = arg 
b 
04 — tan? (-) 
a 
d 
0, = tan”! (=) 
Cc 


Seja: 





al 
c+di 


senh[(1 + )Bx] 
E 

senh[(1 + BL] 

eU+DBx o eTU+DBx eBxpiBx o e Bxe-ibx 
“= CUFDAL- e-CHDPL ePLeiPL — e-PLe-iPL 


* ef*[cos(Bx) + isen(Bx)] — e É*[cos(Bx) — i sen(Bx)] 
“E BL [cos(BL) + i sen(BL)] — e-ÉL[cos(BL) — isen(BL)] 


* cos(Bx) [e?* — e-Bx] + isen(Bx) [eP* + e-bx] 
“ * cos(BL) [ei — e-PL] + isen(AL) [eB! + e BL] 


o cos(bx) senh(Bx) + isen(fx) cosh(Bx) 
Ki cos(BL) senh(BL) + isen(BL) cosh(BL) 


Desta maneira, para determinar-se A, tem-se: 
4= [cos(8x) senh(8x)]2 + |sen(8x) cosh(Bx)]? 
— Alcos(BL) senh(BL)]2 + [sen(BL) cosh(BL)]|? 
o 7 (cos2(Bx) [e?*x — e-Bx]2 + sen2(Bx) [e?* + e-Bx]2) 
E 77 fcos?(BL) [e?! — e-BL]2 4 sen?2(BL) [e?! + e-BL]2) 
ae (cos2(Bx) [e2bx — 2eBxe-Bx + e-26x] + sen?(Bx) [e28x + 2eBxe-Bx + e-26x]) 
— dfcos?(BL) [e2BL — 2eBLe-BL + e-2BL] 4 sen?(BL) [e28L + 2eBLe-BL + e-2BL]) 
= fcos2(Bx) [e28x — 2 + e-26x| 4 sen2(bBx) [e2bx + 2 + e-26x]) 
pus fcos2(6x) [2 cosh(26x) — 2] + sen2(6x) [2 cosh(28x) + 2]) 
— Afcos2(BL) [2 cosh(28L) — 2] + sen2(BL) [2 cosh(2BL) + 2]) 
4= 2 cosh(28x) [cos2(Bx) + sen2(Bx)] — 2[cos?(Bx) — sen?(Bx)] 
— AZ cosh(28L) [cos2(BL) + sen2(BL)] — 2[cos2(BL) — sen2(BL)|] 


Das relações trigonométricas, tem-se que: 
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cos“(y) + sen*“(y) = 1 
cos*(y) — sen?(y) = cos(2y) 
Logo: 


] 
= 2 cosh(26x) — cos(28x) “(2 cosh(26x) — cos(2Bx) /2 
— A2cosh(28L) — cos(28L) |2cosh(28L) — cos(28L) 
Dessa maneira, a fase é dada pela expressão 


o — a) 


ida senh[6L(1 + i)] 


Chega-se, desta maneira, às expressões de A e ly mostradas em Carslaw e Jaeger 
(1959). Portanto, a solução para o regime periódico permanente é: 
6(x,t) = u(x)eitwt+e) = AQo)e"Peilwt+e) 
9(x,t) = AG)etvi++P = AG)[cos(wt + E +'W) + isen(wt + E +'W)] 


Assim, a solução para o problema T(x,t) é dada pela Equação (A.5): 


T(x,t) = A(x) sen(wt + e + !P) (A.5) 


Onde 


— (2 cosh(28x) — cos(2Bx) la 
E E cosh(28L) — DID] 


re senh[Bx(1 + i)] 
E senh[BL(1 + d)] 
Ju 
“N2a 
Neste trabalho, adotou-se: 
da | 
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